The derivative of the topological susceptibility at zero momentum and an 
estimate of rj mass in the chiral limit 



J. Pasupathy"", J. P. Singh^, R. K. Singh"^ and A. Upadhyay^ 

Center for High Energy Physics, Indian Institute of Science, Bangalore-560 012, India 
* Physics Department, Faculty of Science, M.S. University ofBaroda, Vadodara-390 002, India 



in 
o 
o 

(N 

Oh' 
(D 

m 
(N 



> 
O 

\o 

(N 

a^ 
o 
vn 
o 



X 



Abstract 

The anomaly-anomaly correlator is studied using QCD sum rules. Using the 
matrix elements of anomaly between vacuum and pseudoscalars vr, rj and t]', 
the derivative of correlator x'(0) is evaluated and found to be 1.82 x 10~^ 
GeV^. Assuming that x'(0) has no significant dependence on quark masses, 
the mass of -q' in the chiral limit is found to be ss723 MeV. The same calcula- 
tion also yields for the singlet pseudoscalar decay constant in the chiral limit a 
value of PS 178 MeV. 



In QCD, it is known that the vacuum gauge field configurations are significant, in particular the instanton 
solutions corresponding to self-dual fields (GJ^,^ = ±(1/2)6*^;^) play a role in solving the so called U{1) 
problem, i.e., the ninth pseudoscalar rj' remains massive even in the chiral limit when all quark masses 
are zero. The axial vector current in QCD has an anomaly 



d'^qi.l^q = 2im, q^^q - ^ G;l,G'"'^^ where, G'^^^^ = \e>^'""'G%. 
The topological susceptibihty xio^) defined by 



a 



d^x e*^-"(0|r{Q(x), Q(0)} |0), with, Q{x) = ^ G^.G'^'^'^ 



Svr 



(1) 



(2) 



is of considerable theoretical interest and has been studied using a variety of theoretical tools like latice 
guage theory, QCD sum rules, chiral perturbation theory etc. In particular the derivative of the suscepti- 
bihty at = 



x'{0) 



dq' 



(3) 



enters in the discussion of the proton-spin problem ni|2l|3l0|3. In the QCD sumrule approach, one can 
determine x'(0) follows. Using despersion relation one can write 



X'{q') X'{0) 



vr 



ds Q{x{s)) 



+ 



s{s — q^Y ^'^{^ ~ q^) 



+ subtractions . 



(4) 



Defining the Borel transform of a fuction f{q^) by 

Lim 



Bfiq') 
one gets from Eq.@ 

x'(o) 



-q , n — > oo 



('_„2\n+l / rl \ n 



(5) 



X'{q'' 



(6) 



According to Eq.©, '^{x{s)) receives contribution from all states |n) such that (0|Q|n) ^ 0. In partic- 
ular we have L6J 

{0\Q\n') = ^^ml (^^^^) (7) 
\md + muJ 2V2 

The matrix elements, when |n) is |?7) or |ry'), can be determined as follows. It is known from both 
theoretical considerations based on chiral perturbation theory as well as phenomenological analysis that 
one needs two mixing angles and 9q to describe the coupling of the octet and siglet axial vector 
currents to rj and rj' (71|H]|51. Introduce the definition 



(0|j;5|P(p)) =i/^p^; a = 0,8; P = ??,ry', (8) 

fi5 



8 

where J 5 are the octet and singlet axial currents : 



•^^5 = {u^til5U + djf,-f5d - 2s7^75s) (9) 

= -^{uit,l5U + djf,-f5d + SJf,j5s) . (10) 

The \P{p)) represents either rj or rj' with momentum p^. The couplings fp can be equivalently repre- 
sented by two couplings /§, /o and two mixing angles 6*8 and 6*0 by the matrix identity 



V 4' / \ fs sin 08 fo cos 6^0 



(11) 



Phenomenological analysis of the various decays of t] and 77' to determine fp has been carried out by a 
number of authors QIEIIll- In a recent analysis ||9| Escribano and Frere find, with 

/s = 1.28/^ (A = 130.7MeV), (12) 

the other three parameters to be 

08 = (-22.2 ±1.8)°, 00 = (-8.7 ±2.1)°, /o = (1.18 ± 0.04) (13) 
The divergence of the axial currents are given by 



d^'J^ = {ruu wy^u ± dj^d - 2ms s'y^s) (14) 
9^J°5 = ^{m^uj,u + mddj,d + msSy5s) + ^^G%G''f''' (15) 
Since niu, nid << nis, one can neglect them 1 10| to obtain 

{0\-^G%G''''"\v) = ]J-ml [fs cos es-V2fosmeo) (16) 
(0|^G^,G"^^n^') = (/8 sin08±V2/o cos0o). (17) 

Using Eqs.Q, (IT6l and (fTTl we get the representation of xio^) in terms of physical states as 

4 /™ ^ \ 2 ™4 



= -^r^ 2\ \ ~ oAf 2 — 27 [fs ^°^^s ~ ^2/0 sm0o 



- r, — 5— /8 sin 08 ± V 2/0 cos 00 ± higher mass states. (18) 
24(g^ - m^i) V / 



On the other hand, has an operator product expansion IlllfT^ fTlBI 



x{q'^)oPE 



a. 



lb TT TT 



(0| — G^IO)^ + 16 

TT 



128 g4 

,4 




y dp n(/9) K|(Q/9) + screening correction to the direct instantons 



.(19) 



In Eq.(fT9t. the first term arises from the perturbative gluon loop with radiative corrections fT2ll . the 
second, third and fourth terms are from the vacuum expectation values of G^, G'^ and G^. The (0|G^|0) 
term has been expressed as (OlG^lO)^ using factorization 1 11 1. The fifth term proportional to the quark 
mass has been computed by us and is indeed quite small compared to other terms numerically. Finally, 
the last two terms represent the contribution to xiq^) from the direct instantons flTTl . n{p) is the density 



of instanton of size p, K2 is the Mc Donald function and Q 



-q . In a recent work [13], Forkel has 



emphasized the importance of screening correction which almost cancels the direct instanton contribution 
(cf. especially Fig. 8 and Sees. V and VI of Ref. (T3i ). For this reason we shall disregard the direct 
instanton term and screening correction for the present and return to it later. 



From Eq.©, we now obtain 



x'(o) 



nid + mu 



1 + 



M2 



e A/2 



+ TT7 f/s cos 6^8 - \/2/o sin Oq 



24 




e m2 



+^ (/s sin 08 + ^/2/o cos 0o) ' + 



e 



a 



4^) -,M'E,{W'/M^ 



-16 



3 



as 74 Us 9 
1 + — — + — - 
vr 4 vr 2 



+ 



,47r / 

I I as 
16M4 V 



1 ^ /- \ 9 
2^ miqiqi) 



1 



M2 



i=u,d,s 



64 M2 



In 



^^2 



M2 



a 



IT 



vr^ as 



128 M6 vr \ TT 



G' 



(20) 



Here E'o(x) = 1 — e ^ and takes into account the contribution of higher mass states, which has been 
summed using duality to the perturbative term in xope, and W is the effective continuum threshold. We 
take = 2.3 GeV^, and Fig.l plot the r.h.s. of Eq.(|20|l as a function of M^. We take as = 0.5 for 
p = I GeV and 



Writing 



(Ol^^G^lO) = 0.5 GeV^, (0|ss|0) = 0.8(01^^10) with {0\uu\0) 
nis = 150 MeV and mu/md ~ 0.5. 



{0\g!G'\0) = -{0\g'sG'\0), 



-(240 MeV)' 



as in Ref. (5], we take e = 1 GeV^. We also have the PC AC relation, 

-2{mu + md) {0\uu\0) = ml. 



(21) 



(22) 



(23) 



For /o, fs, Og and 6*0 we use the central values given in Eqs.lfT^ and ( fT^ . 
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Fig. 1: Various ters contributing to x'(0)' Eq. <20> . The value of x'(0) is the one obtained without the direct instantons. The 
latter, see Eq.(29), is given by x'di^ which is larger than x'ope and also has the wrong A/^ behaviour suggesting that screening 
corrections are important. 



Let us now examine how the various terms in the r.h.s. of Eq. d20b add up to remain a constant. The 
pion term is small and has little variation because of its low mass, rj and rj' are significantly larger and i] 
is even larger that 77'. In Fig.l the upper Une gives the combined contribution of vr, 77 and 77' which we 
denote as x'poies ^^^^ ^^^'^ ^^^^ ^ gentle increase with M^. The OPE term given by the last three 
lines in Eq. (l20t . which we denote by x'op^;' that 

X'(0) = Xpoles ~ 'X'oPE 

is also plotted in Fig.l. It is seen that x'ope roughly about 25% of x'poies ^^^^ increases with M^, with 
the result that x'(0) is nearly a constant w.r.t M^. 

We expect this trend of compensating variation in x'poies ^'^^ x'ope maintained when variation 
in x'poies to uncertainties in ^g, /s; /o [see Eqs.(IT2ll and ([131 and the variations in x'ope '^^^ 
to uncertainties in the estimates of the vacuum condensates are taken into account. We can then obtain, 
from Fig. 1, the value 

x'(0) w 1.82 X 10"^ GeV^. (24) 

We note that the above determination, Eq. (l24li . is in agreement with an entirely different calculation by 
two of us L14J from the study of the correlator of isoscalar axial vector currents 

n^=° = Ij e^'"(0| {nif.l^u{x) + d^^l^d{x),u^^-i^u{Q) + ^7^75^(0)} |0) 

n^=° = -Il{=\q^)g,,,+Iii=\q^)q^q,. (25) 

I\.{^^{q^ = 0) can be computed from the spectrum of axial vector mesons. In Ref. fT4ll a value 



Yi{-^{q^ = 0) = -0.0152 GeV^ 



(26) 



was obtained. It is not difficult to see that when m„ = = 



x'(o) 



1 



Ui=\q' = 0) 



(27) 



which is consistent with Eqs. d24t . Let us now return to Ea.(fT9l and consider the effect of incorporating 
the direct instanton term in Eq. d20t in the spike approximation |^: 



n{p) = no5{p- pc) 



(28) 



with no = 0.75 x 10 GeV and pc = 1.5 GeV . The contribution of the direct instanton to 
B[y^{q^)/q^] can be found using the asymptotic expansion for K2{z) and i^2(^) ^'^'^ ^^^^ 



9 1 45 
Mpc + -,^^ + 



4 Mpc 32 M^p 



(29) 



We have plottend this term separately in Fig. 1 . We note that unlike Xpoies ^'^'^ x'opE' which increase with 
and therefore compensate each other, the contribution of x'di^ Eq.d^. decreases rapidly with K'P. 
It is not difficult to see that x'(0) will no longer remain constant. This strongly suggests that screening 
corrections to x'di important just as they are for [x{(f')/'f'] as found by Forkel Iil3i . 



We now turn to an estimate of -q' mass in the chiral limit: 



0. In this limit 



5C/(3) flavor symmetry is exact and, we have m^r 



while T]' is a singlet. Let us denote by 



— ^'("T-s = 0) arid = m^i (m^ = 0), the singlet particle and its mass in the chiral limit. Returning 
to Eq.(fT9l. we first note that the exphcitly quark mass dependent term vaxoPE 



■16 ( ^) E Mmi) ~ 1-9 X 10-6 GeV^ 

i=u,d,s 



is numerically much smaller than, for expample 



_9_ /a, 
64 V vr 



4.5 X 10"^ GeV" 



which itself is much smaller than the perturbative term. In the chiral limit (OlQIvr) = {0\Q\r]) = 0. If we 
assume that the quark mass dependence of x'(0) is neglig ible then x'(0) in Ea.(l20l can also be expressed 
in term of /^^ and as: 



x'(o) 



From Eqs. d20b and d30b we may then write 




e 



B 



(30) 



M2 



+ 



-f (1 + ^ 



1 

24 



e 



(/s cos 6*8 - \/2/o sin 




e 



+ — (/8 sin 08 + ^^/0 cos 0o) 1 + ^ 



e M-' 



(31) 



We can find /^^ and from Eq. d3Tl using the least "chi-squared" criterion in the range 0.8 GeV^ < 
Af2 < 1.5 GeV2. We find 723 MeV and corresponding fr^^ = 178 MeV. In Fig. 2 we have plotted 

the l.h.s. and r.h.s. of Eq. dSTt as a function of AP for best-fit values of and /^^. The decay constant 
fr,^ is of the same order as physical decay constants /s and /q. 
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Fig. 2: Estimate of r] mass and coupling in tlie chiral limit, see Eq.(31). Thie continuous curve corresponds to = 723 MeV. 



We now compare our result for x'(0) with some earlier results. In Ref. fll. Narison et al obtained 
a value for x'(0) ~ 0-7 ^ 10~^ GeV^, substantially different from the value derived here. Since the 
expression for xope used by us is identical to theirs, albeit the estimate used for the gluon condensates 
are slighly different, we need to explain the difference in the end result for x'(0)- The most important 
difference is in the expression of x{q^ ) terms of physical intermediate states. We have seen that both -q 
and rj' contribute, and in fact rj makes a larger contribution than t]'. In Ref. 1 1 1 only ?7'(958) state is taken 
into account. We have also seen that if we were to take the chiral limit then rj and rj' contribution to 
is representable by 77^ with mass = 723 MeV, which is substantially different form the physical r/' 
mass. This also explains why Narison et al. find stability in the sum rule only for rather larger W^(=6 
GeV^) instead of our W^=2.3 GeV^. We must also add that while our Eq.© involves only 
Narison et al. use the linear combination of two sum rules (cf. Eq.(6.22) of Ref. 11 1). 

In Ref. [2 1, loffe and Khodzhamiryan's claim that the OPE for xil'^) does not converge is based on 
the following. They computed the correlators 



d'x e^'^-^0\T{j%ix),4,iO)}\0) 



where, J^g = ^7^75^ {q = u, d, s) with m„ 
current. Introducing the definition 



rud 



but rus / and j^^ is the flavor singlet axial 



they estimated 

~ 2.5. (32) 

If SU{3) symmetry were to be exact, this ratio would be unity. Insisting that the ratio in Eq. (l32t should 
be close to unity even when rUs 7^ 0, they concluded that their result signals a breakdown of OPE |2 1. As 
discussed earlier, (0| J^5|77') / 0. In fact using the values given in Eqs. (fT2t and (fT3l . it is easy to obtain 

9?,' \/2(/o coso -V^fs sin 6*8) 



/s sing +\/2/o cos 60 



2.24 



(33) 



which is close enough to the estimate of Ref. fT|. As in the case of Narison et al [1 1, loffe and Sam- 
sonov [51 and Forkel [131 also do not take into account the vr, rj matrix element of the anomaly in their 
sum rules involving x{q^)- 

In conclusion we find a value of x'(0) f« 1.82 x 10^^ GeV^ without incorporating direct instantons. 
Screening corrections to the latter appears to be significant. We also obtain an estimate = 723 MeV 
and = 178 MeV. 
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